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1. Introduction 

Eigenvalue problems with matrix Hamiltonians play essential role in many physical applications. 
For instance, in scattering theory of composite particles they appear as multichannel 
Hamiltonians. In quantum optics they are used for describing multilevel atoms interacting 
with a quantized radiation field. In the simplest case of a two-level atom this is a 2 x 2 matrix 
Hamiltonian and if the rotating wave approximation is used the well-known exactly solvable 
version of the Jaynes-Cummings (JC) Hamiltonian is involved (see e.g. jU I2])- Sukumar 
and Buck presented a number of exactly solvable models of a two-level atom coupled to 
single mode radiation. Here we are addressing the question if this class of exactly solvable 
Hamiltonians can be enlarged with the help of the well-known supersymmetry approach (or in 
other words, by means of Darboux transformations). 

While analyzing this question we were aware that the usual multichannel supersymmetric 
quantum mechanics was almost useless for this purpose. The reason is that this procedure is 
developed under the additional condition that the superpotential matrix W{x) be Hermitian 
|1] but, as we will show below, this is not just the case of the Jaynes-Cummings Hamiltonian. 
In this respect Andrianov et al already considered the case where this assumption could be 
replaced by a weaker demand, namely, the derivative of the matrix W{x) is Hermitian, which 
is equivalent to W{x) — W^{x) = C, where C is a skew-hermit ian constant matrix. However, 
they excluded the possibility for C to be proportional to the Pauli matrix a2 (thus making 
impossible to apply their method to the Jaynes-Cummings Hamiltonian) and moreover, their 
approach was developed in detail only for a particular case. So, the main aim of this paper is 
to analyze carefully the problem of hermiticity of the potential difference produced by general 
Darboux transformations |6ji7j and then apply these results to the Jaynes-Cummings eigenvalue 
problem. 

Another interesting point we would like to stress here is a new character of the factorization 
property that we observe when dealing properly with the hermiticity problem. To show this 
feature, first we find the transformation operator L, defined according to Goncharenko and 
Veselov |6j as an intertwiner of two matrix Hamiltonians: Lho = hiL. Once it is known, 
we can find its formally adjoint L'^ and compose it with L which gives us L~^L = ho + S, 
where S, in contrast to the usual SUSY approach is not proportional to the identity 
matrix but, in general, it is a nontrivial (matrix) operator. Since L"*" takes part in the adjoint 
intertwining relation, L^hi = hoL^, the superposition L~^L is a symmetry operator for ho 
meaning that actually is a symmetry operator too. Such an appearance of symmetry operators 
was first noticed in [Hj and treated as a hidden symmetry exhibiting within supersymmetry. 
Unfortunately, because of the above mentioned restriction imposed on C, the method of [3] 
cannot be applied to get the hidden symmetries for the Jaynes-Cummings Hamiltonian, so that 
in this repect we will make use of our more general set up. 

The paper is organized as follows. In the next section, to fix the notations, we briefly describe 
the solution of the multiphoton Jaynes-Cummings eigenvalue problem pointing out how one 
can find "non-physical" solutions. Section 3 is devoted to adapting the known technique of 
matrix Darboux transformations to Jaynes-Cummings Hamiltonians. In section 4 we give a 
careful analysis of general matrix Darboux transformation leading to Hermitian potentials. In 
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particular, we find that they allow for a new factorization scheme responsible of appearance of 
non-trivial symmetry operators. We indicate a way for establishing a one-to-one correspondence 
between the spaces of solutions, and finally construct an integral transformation operator. 
Section 5 establishes the link between the results of the previous section and the Jaynes- 
Cummings eigenvalue problem. In Section 6 we explicitly construct new exactly solvable matrix 
Hamiltonians of the Jaynes- Cummings type. Some conclusions are drown in the last section. 



2. Jaynes-Cummings eigenvalue problem 



Our starting point is the matrix Hamiltonian 

_ f N + a pa'' \ , . 

describing a multiphoton interaction between a two-level atom and a single mode radiation 
field El- Here a and are the standard harmonic oscillator creation and annihilation 
operators, N = a~^a is the number operator, a G M and /5 G C are parameters. The real integer 
k corresponds to the number of photons the atom exchanges with the field. For the most 
interesting particular case of the one photon exchange, k = 1, we will use the special notation 
hjc '■= h^^^ . In what follows we need not only the known "physical" eigenfunctions of of hjc 
corresponding to its discrete spectrum eigenvalues but also those solutions of the Schrodinger 
equation which do not belong to the Hilbert space. Therefore, below we give a short outline of 
how these solutions may be obtained. 

Let ipn be the usual number operator eigenstates, Nipn = nipn, with the properties 

i" 1/2 _l_ i* ^ / '2i 

a ipn = + 1 ~ ^)fe ^n-fc, {(^'^) Wn = + "^n+fc- Here {x)k denotes the standard 
Pochhammer symbol, {x)k = r(a; + k)/T{x). The Hamiltonian ^ can be diagonalized by 
a rotation inside any two-dimensional space spanned by the vectors (V'njO)* and {Q,ipn+kY- 
(The superscript "t" meaning the transposition.) Indeed, if we write an eigenvector ^! e oi h^''\ 
h^'^^'^E = E^i El in the form e = {ciipn, C2'?/'n+fc)*, then the eigenvalue problem is reduced to a 
simple system of two linear homogeneous equations for the coefficients ci and C2, which gives 
us the spectrum 



E = E'/ = n + l{k + a)± ^\{a - kf + \P\\n + 1), (2) 
and the coefficients 

C2 = TTT , Ci = \l\ — ICoP . (3) 

/3 [|/3P(n + l)fe + (n + a-E)2]i/2 

The set of eigenvectors thus determined, \E'„ = {cx'^n-, C2V^n+fc)*, 'n = 0, 1, . . . is not complete 
in the Hilbert space but the missing eigenvectors can readily be obtained by noticing that 
the vectors \E'^ = {0,ipjy are eigenvectors of h^''\ corresponding to an equidistant part of the 
spectrum h^'^^'^j = j^j, j = 0, . . . , k — 1. 

An interesting feature of the Hamiltonian is that for k > 2 the square root in (j21) 
dominates for large n leading to the unbounded character of the spectrum both from above 
and from below. From this point of view this Hamiltonian is similar to a Dirac Hamiltonian. 
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The easiest way to get "unphysical" solutions for hjc is to notice that in the above 
construction we needed just a couple of states V'n+fc} related by the operators a'^, {a^Y ■ 
We can take also these vectors from the unphysical sector. 

For instance, choosing them among the basis vectors {ip-n}'^=i of the skew-hermit ian 
representation of the oscillator algebra where a'^'ip-n = a/— n + 1 V'-n+i, a'tp-n = \/—nip_n-i^ 
a'^ip-i = we will get eigensolutions \E' = (ci'?/'_„, C2'>p-n+k) corresponding to an eigenvalue 
E. The expressions for E and ci,C2 are given in this case by the same formulae (j21) and 
(121), where n has to be simply replaced by —n. In the coordinate representation where 
ipn{x) = (A/7rn!2"')~^/^ exp(— a;^/2)if„(x), Hn{x) being the Hermite polynomials, the negative 
eigenvectors are realized as (x) = ipn-i{ix). Here the eigenvalues can be complex according 
to the sign under the square root in Q that now may be negative depending on n, k, a, and 
(3. A particular set of nonphysical eigenvectors is given by [ip-j, 0)*, j = 1, . . . , k — 1. 

Even we can use a general solution C,i{x) of the oscillator equation a~^aC,i{x) = eC,i{x) 
together with ^2 = (ct"^)^^!, provided ^1 oc a'^^25 so that the nonphysical eigenfunction has the 
form E = (ci^i, 02^2)- The coefficients ci, C2 depend now on e and can be easily computed, but, 
in general, they are complex too. 



3. Darboux transformation of the Jaynes-Cummings Hamiltonian 



To construct the Darboux transformation for the Jaynes-Cummings Hamiltonian we are using 
the existing approach developed for the matrix Schrodinger equation |BJ[3. Therefore, we first 
reduce the Jaynes-Cummings eigenvalue problem to the Schrodinger equation with a matrix- 
valued potential. Then applying Darboux transformations we get new exactly solvable matrix 
potentials for the Schrodinger Hamiltonian and the use of the inverse transformation allows us 
to obtain desired transformation for the Jaynes-Cummings Hamiltonian. 

We start with rewriting hjc in the coordinate representation as 

ho = -dl + Vo + b-fd, , 1= J ^ , hjc = \ho (4) 

with 

where, for simplicity, we have taken (3 real. The operator contains the first order derivative 
term, therefore the method of [HI IZ| cannot be applied to it directly. This undesirable term is 
easily removed by the unitary transformation 

JiQ = U+/ioU , U = (U+)^i = e^^"" (6) 

which gives us the Schrodinger Hamiltonian 

ho = -dl + Vo, Vo = U+VoV~\b\ (7) 
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The Darboux transformation operator L for a Schrodinger type Hamiltonian is known to be 
in the form j^l [7] 

L = d,-W, W = uJa-\ (8) 

leading to the transformed potential 

V^ = V,~2W^. (9) 

Here and in the following the subscript x means the derivative with respect to x and the 
tilde over a symbol marks a quantity related with the Schrodinger equation, unless otherwise 
indicated. The matrix W is a solution to the matrix eigenvalue problem 

hoU = Uh (10) 

with eigenvalue L being a diagonal matrix, L = diag(li, I2). In this case Ik are eigenvalues of 
corresponding to column- vectors Uk = (wi.fc, ^2,^)*, = 1, 2, and the matrix U is just composed 
of these columns, U = {Ui, U2)- 

Once the Hamiltonian hi = —d^ + Vi is determined we realize the inverse transform to 
express all quantities in terms of the eigenvalue problem for h^. For instance, U = IJU will be 
an Hq matrix eigensolution, preserving the same eigenvalue A, 

hoU=UL. (11) 
As for the potential Vi, defined by the expression hi = U/iiU^ := —d"^ + Vi, we get 

Vi = Vo- 2W^ + b{-fW - W-f) . (12) 
It is easy to see that the operator L (jH)) is covariant under the transformation U: 

L = UZU+ = d^-W. (13) 

Keeping the usual SUSY terminology (see e.g. jl]) we will refer to the matrix W := UJA~^ 
as a Jaynes- Cummings superpotential. It is interesting to notice that the change in the 
transformed potential (fT^ has two contributions. The term —2Wx in (fT^ corresponds to 
the usual expression for the potential difference due to the Darboux transformation (see e.g. 
[HI El)- The next term proceeds from the first derivative member of the Hamiltonian ho and 
coincides with the potential difference for Darboux transformed Dirac systems (see ^^). 

We see from ()12|1 and (fT!?|) that the new potential and the transformation operator are defined 
in terms of a matrix-valued solution, U, of the initial Hamiltonian. Therefore, as usual we 
call it transformation function. To keep the potential difference regular we have to impose an 
additional condition on U: detU 7^ 0. Another additional condition must be imposed on U if we 
require that Vi be represented by a self-adjoint matrix. It is not difficult to see that if we restrict 
ourselves to real potentials, the last term in (fT^ is symmetric for any matrix W. Therefore, the 
potential Vi will be self-adjoint if W is real and its derivative symmetric, i.e., W — = Cjc 
where Cjc = ~C~^c ^ skew Hermitian constant which because of the real character of W 
should be proportional to 7, Cjc = C07. Now using the relation W = UWU'^ + I7 we find 

W-W+ = c-f, c = Co-b. (14) 
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We notice that both superpotentials W and W can never be simultaneously Hermitian. In 
sectionElwe illustrate that both possibilities (the hermiticity of W or that of W) may take place. 
From here it follows that the usual approach of the multichannel supersymmetric quantum 
mechanics based on Hermitian superpotential matrices, W = W^, cannot by applied to the 
Hamiltonian The above condition should be replaced by a weaker demand ()14j) . 

So, we have established how the matrix Darboux transformation should be modified when 
applied to the Jaynes- Cummings Hamiltonian. It can be considered as a special case of the 
general transformation introduced by Goncharenko and Veselov jH] with an additional condition 
of type ()14p. Now we are going to investigate implications this restriction imposes in general and 
then we will apply obtained results to the particular case of the Jaynes-Cummings Hamiltonian. 

4. Properties of the Darboux transformation for the matrix Schrodinger equation 

To make easier the notation in this section we will suppress all tildes over the quantities related 
with the Schrodinger equation introduced in the previous section. 

Darboux transformation of the usual (one-component) Schrodinger equation has a number of 
nice properties (see e.g. [H]) making it popular in different fields of theoretical and mathematical 
physics. In this section we show how some of them are translated to the matrix level. In our 
opinion the most interesting properties are the following: 

• The most convenient way to introduce the transformation operator associated with the 
name Darboux is to define it as a differential intertwiner between two Hamiltonians ho and 
hi. Defining (formally) an adjoint operation such that the Hamiltonians are self-adjoint 
one gets the adjoint intertwining relation 

Lho = hiL , L+hi = hoL+ . (15) 

• Given L and L"*" one can compose them thus discovering factorization properties first 
discussed by Schrodinger and studied in detail by Infeld and Hull. 

• In matrix notation the above properties can equivalently be rewritten as commutation and 
anticommutation relations giving rise to a supersymmetry algebra commonly known also as 
supersymmetric quantum mechanics. 

• There exists a procedure allowing one to realize a one-to-one correspondence between the 
spaces of solutions of the Hamiltonians /iq and hi. Therefore, if /iq is exactly solvable, 
then hi is solvable as well. In this case it is rather easy to find the spectrum of hi when 
the operators are defined in a Hilbert space. In particular, one can conclude that the 
spectrum of hi coincides with the spectrum of h^ with the possible exception of a finite 
number of levels. In the simplest case either the spectra may differ only by the ground 
state level (so called case of the "exact supersymmetry") or they may be identical ("broken 
supersymmetry" ) . 

• Making use of a special two-step differential transformation one recovers an integral 
transformation (sometimes called Luban-Pursey method) leading, in general, to a relation 
between the method of Darboux transformations and inverse scattering method |12i| . 
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• There exist nice Crum-Krein formulae for a compact representation of the resulting action 
of a chain of first order transformations. Recently they have been generalized to the matrix 
case so we will not discuss this point here. 

4.I. Factorization properties and underlying supersymmetry 

First we will study how the factorization properties are translated to the matrix level. 
Since 

W+ = W-C, W =UJA'^ (16) 

hdA = Uh (17) 
we have = —d^ — W + C which results in the following superposition 

L+L = ho + go, go = CL- UhU-^ . (18) 
To find the superposition LL^ we express L in terms of W^: L = — — C which yields 

LL+ = hi+gi, gi = - {UUA-^Y . (19) 

From here it follows that g^ = g^ commutes with and gi = g^ commutes with hi and they 
are intertwined by L, Lg^ = giL. The facts which can also be checked by the direct calculation. 
Nevertheless, it is necessary to stress here that this observation is correct only formally. When 
all operators are considered as acting in a Hilbert space they may be not commuting because 
of different domains of definitions and the subject should be studied more carefully. 

From ()1U|) and its adjoint form it follows another nice property of the transformation function 

U: 

2{UUA-^)^ = CVi - VqC = Chi - hoC 

which shows that for the particular case C = the matrix WLW^^ is constant. We would like also 
to stress that even in this particular case the factorizations fll8|) and (fT^ do not coincide with the 
ones giving rise to the usual multichannel supersymmetric quantum mechanics: L^L = h^ — XI 
and LL^ = hi — XI where A is known as factorization constant and / is the unit matrix. 
This fact was already mentioned in |7|. Another point worth to be noticed is that in contrast 
with the paper where similar symmetry operators have been found, we give their explicit 
expression in terms of a solution U of the initial eigenvalue problem. In fact such nontrivial 
symmetries may exist only if the Hamiltonian ho has at least one matrix solution U with the 
property (fT^ . 

Another interesting observation is that if the spectrum of ho is nondegenerate and essentially 
self adjoint operators ho and go have a common set (dense in the Hilbert space) in their domains 
of definition, where they commute, they should be related by a functional dependence. This 
statement follows from the fact that if the spectrum of ho is nondegenerate then a complete 
set of operators commuting with ho consists of only one operator which is ho itself. Hence, by 
definition of a complete set of operators (see, e.g [TT]) any other self adjoint operator commuting 
with ho is a function of ho- In this case all eigenfunctions of ho can be found by solving the 
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eigenvalue problem for gQ, which is a first order differential operator. We shall see latter that 
just the Jaynes-Cummings Hamiltonian presents a nontrivial example of such a situation. 

From the point of view of supersymmetric quantum mechanics we can construct here a wider 
(extended) superalgebra than the one usually appearing in the scalar case. In addition to the 
superhamiltonian 



Now, as usual, the intertwining relations ()15p are equivalent to the commutation of H and G 
with the supercharges and the factorization properties are translated to the anticommutation 
relation for the supercharges QQ~^ + Q~^Q = H + G. It is also evident that H and G commute 
between them. According to the remark made in the previous paragraph, these relations are 
also rather formal and when a concrete Hilbert space is considered it is necessary to take care 
of domains where the operators act. 

In the case of 2 x 2 uncoupled Schrodinger equations with Vq = the conjugation property 
for the superpotential can readily be analyzed. We have found that the condition (fTT)|) can be 
satisfied only for C = meaning that differential first order symmetries do not exist in this 
case. On the other hand, the result of the previous section shows that matrix Hamiltonians 
with such symmetry operators do exist. Another useful remark concerns the usual multichannel 
supersymmetry approach. It imposes the restriction on L to be diagonal with equal elements. 
In our case L is allowed to have different diagonal elements which enlarges considerably the set 
of exactly solvable partners for a given Vq. 

4-2. One-to-one correspondence between the spaces of solutions 

Being differential, the operators L and L"*" have nontrival kernels. Therefore, the correspondence 
between the spaces of solutions of the Schrodinger equations for ho and hi should be studied 
carefully. More precisely, a detailed analysis is necessary for kernel spaces of these operators 
which are eigenfunctions of ho and hi corresponding to the eigenvalue matrix L. In the scalar 
case Keri^ho — /) is a two dimensional space easily determined by a solution u, hou = lu. This 
possibility is based on the property of the Wronskian of two solutions corresponding to the 
same eigenvalue to be equal to one. First we are going to show how this property translates to 
the matrix case and then how it helps us to establish a one-to one correspondence. 

The conjugation equation for the superpotential being rewritten in another form 




and mutually conjugated nilpotent supercharges 




we have the symmetry operator 




W{U, U) := U+U - U+U^ + U+CU = 



(20) 
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establishes a property of the transformation function U. We relate it with the fact that the 
Wronskian of a solution with itself vanishes. Therefore we call the quantity W{U,U), defined 
by the middle part of ()20|) . Wronskian with coinciding arguments. Quite naturally for different 
arguments it is defined as: 

W{V,U) ■.= V+U -V+U^ + V+CU . (21) 

We notice that for C = and U, V being usual functions (not matrices) it coincides with the 
usual expression for the Wronskian. Therefore, we consider fl21|) as a matrix generalization of 
Wronskian. 

Now we will show that if W is a solution to (jl7|) . satisfying then there exists another 
solution to the same equation, V, which can be found from the condition W{V,U) = —1. 

We proceed first to find an explicit expression for V. Multiplying the equation W^(V, U) = —1 
hj from the right and using we get 

V+-V+(W+)-iW+=W-^ (22) 

Multiplying the adjoint form of ((22) by U'^ from the left we obtain (U'^V)^ = W"^(W+)~^ 
which gives us the final result 



V = U I {U+U)-^dx + UCi. 



(23) 



In general, we cannot neglect the matrix integration constant Ci. This is due to the evident 
property of ()17|) not to be a linear equation with respect to the multiplication on matrices. 
Only if this constant takes a definite value this function is a solution to (fT7|) . 

Now we proceed to find Ci. For this purpose we put the function into equation (fT7|). 
Using the same equation for U and (fTBj) we obtain 

PX PX 

VoUCi - Ci{U^)-^ - U^^Ci = u {U^Uy^dxL -uL {U+Uy^dx + UCih . (24) 

Jxo Jxo 

Here it is necessary to calculate the commutator of L with the integral at the right hand side 
of ()24|) but first we have to find the commutator of L with {U^U)~^. 

From (fT^ and (fT7|) it is easy to get 

UA+U - U+UL = {U+CU)^ 
which can readily be transformed to 

which gives us 

PX PX 

/ {u^uy^dx L-L (u+uy^dx = -u-^c{u+y^ + Ca . 

Jxo Jxo 

Here 

C2 = {U-^C{U^y%=,, = -C^ . (25) 
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Now once again using (fTTj) we obtain from the equation for Ci. 

C2 + CiL - LCi = , C+ = Ci . (26) 

Since we just reduced the Schrodinger equation for V to the last equation for Ci, this means that 
with the constant Ci thus determined the function V (j2SI) satisfies the Schrodinger equation 

(HZD. 

Here we would like to notice that for any n-dimensional vector space Ker {Jiq — E) of solutions 
of the eigenvalue equation for with a given E 1^, k = 1,2, . . . ,n, h^tpE = EipE, (we remind 
that /fc correspond to the eigenvectors Uk of the matrix U) the kernel of the operator L is the 
empty set whereas LU = 0. To find solutions of the Schrodinger equation with the eigenvalues 
Ik we can act by L on the function V which yields 

U = LV={U+)-\ hiU = Uh. (27) 

Since L is supposed to be diagonal, columns of U are eigenvectors of hi with eigenvalues Ik- 

An obvious but necessary remark is that L^U = 0. Moreover, it is easy to see that equation 
(fTH|) is covariant under the transformation U hi = (U^)~^ meaning that U has the same 
property. Hence, another solution of the eigenvalue problem for hi corresponding to the matrix 
eigenvalue L can be found using the same formula (j23j) applied this time to U: 

V={U+)-^(l U+Udx + Ci] , hiV = VL. (28) 



Xo 



Once again because of the diagonal character of L columns of V are eigenvectors of hi with 
eigenvalues Ik- The constant Ci from should be determined in a similar way that Ci from 

(El- 

So, for the eigenvectors ipE of hi, hiipE = EipE, with E ^ Ik, k = 1, . . . ,n the kernel of is 
the empty set, KerL"*" = 0. This means that if ipE is also an eigensolution to L~^LipE = ^J'E'^^e 
then (fiE = LipE satisfies the equation LL^LfE = f^EfE- In this case L and L"*" realize a one- 
to-one correspondence between the spaces of solutions foi all E Ik, k = 1, . . . ,n. For E = Ik 
this correspondence may be established by hand U ^ lA, V ^ V and continued by linearity 
at the level of vector- valued eigenfunctions ipE and of /iq and hi. If L+L get ifjE beyond 
the one-dimensional space Span-?/)^, L'^LifjE ^ Span-j/^^;, this correspondence is more subtle and 
needs an additional analysis. 



4-.3. Integral transformations 

Let a matrix constant Cqi be such that (U^)^^Cqi be the solution to the equation (j2Z|)- Then 
the function Vi = V + {U^)^^Cqi with V as given in is the solution to the same equation 
and it can be taken as the transformation function for the next transformation step. After 
some algebra one gets 

{Vi)^iyi)-^ = -W+ + W2, W2:=U 



U+Udx + Co 



(29) 



xo 
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where Cq = Cqi + Ci. So, the potential V2 = Vi — 2(Vi)^(Vi) ^ is given by 

V2 = Vo + AV2 , AV2 = -2{W2). . (30) 

The potential difference AV2 from (jHn|l and is self-adjoint for any matrix U. Nevertheless, 
we can use here only these Us which give rise to a self-adjoint derivative of the superpotential, 
= Wx, since only under the condition (fT^ the function is a solution to the Schrodinger 
equation after the first transformation step. An advantage of this formula with respect to Q 
could be a much easier possibility to get an everywhere nonsingular resulting potential since 
just the constant Cq can be used for this purpose whereas in Q we have not any freedom of 
this kind. Moreover, in such a way one can get families of isospectral and isophase (known also 
as phase-equivalent) potentials. Applying twice appropriately changed formula (jHl) we express 
solutions of the Schrodinger equation with the potential (j3U|) in terms of solutions of the initial 
eigenvalue problem 

ipE = UUA-^iIje -iJEE-{C + W2){^E. - UJA-^iIje) (31) 

h2ipE = <^eE , h2 = -dl + V2 . 

Here may be both vector- valued and matrix- valued eigenfunction of with E being a number 
in the first case and a diagonal matrix in the second case, hoipE = "^PeE. As usual, for E = L 
and ipE = U, the right hand side of vanishes, but the missing eigen-solutions correspond 
to (V+)-i. 

5. SUSY partners for the Jaynes-Cummings Hamiltonian 

5.1. Properties of the Darboux transformation for the Jaynes-Cummings Hamiltonian 

Using the relation between the Jaynes-Cummings Hamiltonian and the Schrodinger 
Hamiltonian given by equation © we can readily apply the results of the previous section 
to the Jaynes-Cummings eigenvalue problem. From now on we restore notations of the section 
El and will distinguish all quantities related with the matrix Schrodinger equation by putting 
tilde over. The quantities without the tilde will be used for the Jaynes-Cummings problem. 

First we notice that according to (fT^ even for a self-adjoint Jaynes-Cummings superpotential 
the constant C introduced in the previous section is not zero but C = by. Factorization 
properties (jl8p and (jl9|) together with the symmetry operators are covariant under this 
transformation. Therefore the symmetry operator for the Jaynes-Cummings Hamiltonian can 
be found from a diagonal superpotential matrix. Since U commutes with C = by the formula 
(j25j) for the constant C2 is covariant also. The unitary transformation (jH)) does not affect the 
eigenvalue matrix L. Therefore the equation for the constant Ci remains intact. The 
equations (P7|) and (PH)) for solutions of the transformed equations with the eigenvalue L are 
also covariant but now they are valid for the self-adjoint Jaynes-Cummings superpotential W. 
For the integral transformation the transformed potential is given by the same formula ()12|). 
where W is replaced by W2, which is given by ()29p. where all quantities are related now with the 
Jaynes-Cummings system. The formula (j3H) for solutions of the equation with this potential is 
also covariant. 
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5.2. Examples 

1. In the simplest case the transformation function can be taken in the form (see section |2I) 

"=(,..% % ). (32) 

It corresponds to the diagonal eigenvalue matrix L = diag(0, 2a— 2) and diagonal superpotential 
W = diag(x, —x). Formula (fT^ gives us the transformed potential 

.r _ ( + 2a - 3 -hx \ , . 

-bx x' + l ) 

and finally the transformed Hamiltonian 

/^, = 2f ^^^ + "-2 -Pa^Y (34) 
\^ —pa a^a J ^ ' 

Solutions of the Schrodinger equation with this Hamiltonian can be obtained from the solutions 
of the Jaynes- Cummings eigenvalue problem using the transformation operator 

(35) 

Using (fH^ we find the symmetry operator 

The spectrum of the Jaynes-Cummings Hamiltonian is nondegenerate and the set of all finite 
linear combinations of its eigenfunctions is dense in the corresponding Hilbert space. The 
operator (yfQo (jH^ is well-defined on this set where it commutes with hjc- Hence, as expected 
the latter is a second order polynomial of g^^. b'^ho = goo{goo -|- &^ — 2a + 2). 

2. Another possibility is to construct the transformation function U from two vectors of the 
type {ciilJn,C2ipn+i) with the coefficients ci,2 given in Q corresponding to the energies E^'"^ 
at /c = 1. For n = we easily get the following superpotential 

i!,). (37) 

We notice that it is singular at the origin and cannot give rise to a self-adjoint potential for the 
spectral problem on the whole real line. Nevertheless, up to a constant shifting it produces the 
same symmetry operator ()36|). Moreover, since at any n both even and odd Hermit polynomials 
are involved in this construction, the superpotentials we can obtain in this way are singular at 
the origin for any n. 

3. As it usually happens in the method of Darboux transformations (see e.g. pi) a 
singular point can be removed by the next transformation if it is realized with the help 
of a transformation function corresponding to the adjacent spectral point. Moreover, two 
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consecutive first order transformations can be replaced by a single second order transformation 
(for details see [7j). As a result we get the following potential difference 



AV 



1 + 2x2)2 + 2x2) 4x^ + 8x2-1) 



We notice here that the spectrum of the Hamiltonian /12 = Hq + AV coincides with the spectrum 



1 2 

of ho except for four levels Eq\ which are removed by this transformations. It is also clear that 



1 2 

taking other values of n one can remove any four levels E^'^_^_i. For this purpose we are using 
two transformation functions 



Wi = 



n{x) \ JJ _ f "^71+1(3^) ■A2lpn+l{x) 
y4i^„+l(x) i)n+l{x) y ' ^ ~ V ^2V^n+2(x) ^„+2(a;) 



and the algorithm developed in [7 . After some algebra we find that the second order super- 
potential is diagonal, W2 = diag(wo,Wi), with the entries Wk = 2ipn+k'ipn+k+i/ {'ip'n+k'^n+k+i - 
ipn+ki^n+k+i) y k = 0,1. Here in the denominator we see the Wronskian of two consecutive 
discrete spectrum eigenfunctions of the usual Harmonic oscillator potential which is known to 
be nodeless [H]. Therefore it gives "good" potential differences for any non-negative integer n. 
For instance, n = 1 case corresponds to the following potential difference 



AV = A 



Ax'^+Sx'^-l L / 

(1+2x2)2 " ys+ix* 1+2x2^ 

4x2 ^\ 8x2(4x-*-9) 

3+4x4 1+2x2 I -L -|- (3+43:4)2 



4. The creation of new energy levels is also possible. For this purpose we need solutions 
of the Schrodinger equation, which do not belong to the Hilbert space discussed in section |21 
They are constructed with the help of the functions ipn{x) = tpn{ix). For the transformation 
function 

U={ ^° 1, n = 0,l, 



!f 

the superpotential reads 

( ^0(3:^) ^^ yo(3:)y>^+i(a:)-yo(a:)yn+i(x) 
ipo{x) ipoix)ip„{x) 
^ 

Using the recurrence relation for the Hermit polynomial we find that its non-zero off-diagonal 
element does not depend on x thus giving us a nontrivial example of a non-Hermitian up to a 
constant superpotential. Substituting the logarithmic derivatives by their expressions we get 
its final form 



W 



X Ay/2n + 2 

PI _|_ 2niHn-i{ix) 
^ H„{ix) 



We notice that for even n the Hermit polynomials Hn{ix) are nodeless whereas if„_i(ix) are 
purely imaginary. So, for n = 0, 2, . . . this formula gives us "good" potential differences provided 
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A is real which is not always the case. For instance, for n = the potential difference is a 
constant diagonal matrix AV = diag(-Bo, —4 — Bq), where Bo = —1 — a ± ^/ {a — 1)^ — 26^. 
Two possible signs here and below are related with different signs in The first nontrivial 
case corresponds to n = 2: 



Bi 



4bx 



l+2x' 

Abx A 4a:-^+3 



—A 4a:-+3 n 



1+2x2 ^ (1+2x2) 



where Bi = — 1 — a ± {a — lY — 66^. We have to notice that after being shifted by a constant 
diagonal matrix this potential difference may be reduced to the previous example at n = 0. 
Hence, to get an essentially new potential difference we have to consider n = 4: 



AV 



rj 86x(3+2x2) 
2 4x4+12x2+3 
8bx(3+2x2) . 16x^+64x6+120x"t+45 



A ibx"+b4x"+i:iUx'+4& D 

— 4 ,.A , , nN9 Do 



4x4+12x2+3 (4x4+12x2+3)2 

where ^2 = -f - a ± ^J{a - If - W. 

5. As the final example we give a one-parameter family of potentials isospectral with Jaynes- 
Cummings Hamiltonian obtained with the help of the integral transformation ()30p . ()29|): 



AV = % 



l+xF(x) L 2x2 , 

F2(x) 4VF(x) 2x-F(x)' 

1_ 2x2 o^ 2x+(x2-l)F(x) 



4'^F(x) 2x-F(x)^ (2x-F(x))2 

where = ^/^[e^^ {c + erf(x)) and real c is such that |c| > 1. 



6. Conclusion 



A careful analysis of the matrix Darboux transformation method has permitted us to establish 
such properties as: (z) a new factorization scheme, which is responsible on appearance of an 
extended supersymmetry underlying matrix Hamiltonians and hidden symmetry operators, (m) 
a one-to-one correspondence between the spaces of solutions, which permits to readily determine 
the changes in the spectrum, and {in) the integral transformation formula, which gives families 
of isospectral Hamiltonians. We applied these results to get a new symmetry operator for the 
Jaynes-Cummings Hamiltonian and build up many of its exactly solvable partners. A future 
step of this research will be the investigation of physical phenomena, such as for instance the 
collapse and revival for new Hamiltonians. 
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